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Abstract 

In classical kinetic or kinetic-like models a particle free path distribution is exponensial, 
but this is more likely to be an exception than a rule. In this paper we derive a linear 
Boltzmann-like equation for a general free path distribution in the framework of Alt’s model 
J. Math. Biol. 9, 147 (1980). In the special case that the free path distribution has at 
least first and second finite moments we construct an asymptotic solution of the equation for 
small mean free paths. The asymptotic solution becomes a diffusion approximation to the 
one-speed Boltzmann-like equation. 


1 Introduction 

Classical neutron transport or transport of thermal energy by photons is described by the 
linear Boltzmann kinetic equation [1-3], where the free path distribution is exponential, i.e., 
collision events form the Poisson process [4]. This is a very restrictive assumption, and, therefore, 
in many cases it is violated [5-8]. Kinetic and kinetic-like models are also used for description 
of transport processes in biology [9-17]. The common assumption in these models is also the 
exponential distribution of the free path, but this is more likely to be an exception than a rule. 

A model, which extends the linear Boltzmann equation to a general free path distribution, 
was proposed in Ref. [18]. This model was embedded in a biological context. The model does not 
take into account absorption and sources, but their inclusion is straightforward. In Ref. [18] a 
partial integro-differential kinetic equation and its diffusion approximation (Patlak-Keller-Segel 
equation) were derived under the assumption that the free path distribution was asymptotically 
exponential. The steady-state case of Alt’s model [18] in a physical context was rederived in 
Ref. [7]. A more general model, including acceleration of the particles due to external forces, 
was proposed in Ref. [19], see also Ref. [20], where a Boltzmann-like equation without absorption 
and sources was obtained. 

Alt’s model [18] can be considered in the framework of the velocity jump model of the 
continuous-time random walks (CTRWs) [9, 21]. A model similar to Alt’s one in the framework 
of the position jump model of the CTRWs was proposed in Ref. [22], see also Ref. [23]. 

In this paper we derive a linear Boltzmann-like equation for a general free path distribution 
in the framework of Alt’s model. In the special case that the free path distribution has at least 
first and second finite moments we construct an asymptotic solution of the equation for small 
mean free paths like in Ref. [24], where an asymptotic solution was constructed for the ordinary 
linear Boltzmann equation. The asymptotic solution becomes a diffusion approximation to the 
one-speed Boltzmann-like equation. In a steady state this diffusion approximation coincides with 
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that obtained in Ref. [25] and differs from the diffusion approximation in Ref. [7] by small terms 
of the second order. 

All the reasoning in this paper is performed in the three-dimensional space. A generalization 
to the n-dimensional space is straightforward. 

The paper is organized as follows. Section 2 describes the model of non-classical particle 
transport. In Section 3 we derive the Boltzmann-like equation, which describes transport with 
an arbitrary distribution of the free path. In Section 4 we construct an asymptotic solution of the 
equation for small mean free paths. In Section 5 the diffusion approximation to the one-speed 
Boltzmann-like equation, as a particular case of the asymptotic solution, is obtained. 


2 Non-classical particle transport: a model 


The model generalizes that of classical linear transport [1, 3], where the extinction and 
scattering coefficients (total and scattering cross sections) do not depend on a free path travelled 
by a particle. In the non-cIassical transport the coefficients depend on the free run path. In this 
case the phase space density of particles takes the form ^(r, v, t, /), where r is a position, r € M^, 
is a velocity vector, t7GVcM^,t>0is time, I >0 is the free path of the particle. The phase 
space density obeys the equation [18] 

d^-^ + vdiC + v-V^ + vaC = 0, ( 2 . 1 ) 


where V means the gradient with respect to r, a = a{r,v,l) is the extinction coefficient (or 
total cross section), v = |r7|. This equation describes noninteracting particles moving at the 
point r with the velocity v so that they stop their free run (scatter or disappear) with the rate 
va. We do not consider here particle acceleration due to external forces; taking it into account 
is straightforward. 

Eq. (2.1) is supplemented by the condition, describing the density of the particles at the 
beginning of their free run. 


i\i=Q^TliT,v,t) = K: 


■ poo 

Jo 


+ F{r,v,t) 


where 


‘^= / K{r,v,v')f{r,v',t)dv' 

Jy 

is the scattering operator, K(r,v,v') = K(r,v' —)• v) is the scattering kernel such that 


( 2 . 2 ) 

(2.3) 


K > 0 and 


K{r, V, v') di) = 1 


[i.e., K{r, -^v') is a probability density function], cjg = ag{r,v,l) is the scattering coefficient (or 
scattering cross section), cr^ < a, F is the source term. The condition (2.2) is a straightforward 
modification of the corresponding condition given in Ref. [18]. This condition means that the 
particles at the beginning of their free run arise as a result of scattering and are generated 
by sources. This condition means also that, if dg < a, a part of the particles disappears (is 
absorbed). 

In general case the phase space density is subject to the initial condition 


^lt=o = 

where iJ is an initial distribution of the density. However, here we impose the initial condition 

^\t=o = X{r,v)6{l), (2.4) 

where 5{-) is the Dirac delta function. This condition means that at the moment t = 0 all the 
particles begin their free run. 
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3 The Boltzmann-like equation 

We define the phase space density 

poo 

4’ir,v,t) = / C{r,v,t,l)dl, 

Jo 

which does not take into account the free path. Eqs. (2.1) and (2.2) yield 
+ V ■ 'V'lp 


(3.1) 


= /C 


poo 1 poo 

’ -v a{r,v,l)^{r,v,t,l)dl + F{r,v,t) (3.2) 

Jo J Jo 


(we assume that ^ ^ 0 as / —oo). 

To express the right-hand side of Eq. (3.2) through V’ we integrate Eq. (2.1) along the char¬ 
acteristics. This yields 


^{r,v,t,l) = 


S{r,v,l)r]^ — , I < vt, 

S{r,v,l)x{'i' — fil,v)6{l — vt), vt<l, 


where 


(3.3) 


(3.4) 


S{r, V, 1) = exp 1“ y a{r — 17(/ — l'),v, I') dl' 
is the survival probability (see Appendix A for details). Therefore, 

poo 

/ cr{r,v,l)^{r,v,t,l)dl 
Jo 

= J p{r,v,l)ri^ — dl + p{r,v,vt)x{r — vt,v), (3.5) 

where 

p{r,v,l) = a{r,v,l)S{r,v,l) (3.6) 

is the probability density function (PDF) of a free path. Eq. (3.5), the Laplace transform and 
the equality h{r — fit) = h{r) yield 


C 


poo 

/ cr{r,v,l)i{r,v,t,l)d 

pvt 

/ P{r,v,l) 
Jo 




I 


r][ r,v,t — - ] dl +p{r,v,vt)e '“*'*^'^x(r,t;) 


^ ^ V, s) + x{r, i’)] 


where 


poo poo 

e"**/(i)dt or Cg{\) = j g{l)dl 

Jo Jo 


is the Laplace transform (we distinguish functions of t and 1). 

Eqs. (3.1) and (3.3) yield 

'ilj{r,v,t) = J S{r,v,l)g^ — f2l,v,t —dl + S{r,v,vt)x{r — vt,v), 
or, equivalently, 

C'il){r,v,s) = ^Cs(^,v, ^ -h i7 • [Cr]{r,v,s) + xir^v)] . 
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(3.7) 


Finally, if we define the memory kernel (f){r,v,l) by [26] 


def Cp{r,v,X) XCp{r,v,X) 
C(j){r,v,X) = - 


which is equivalent to 


we obtain 


CS{r,v,X) l — Cp{r,v,X)^ 

p{r,v,l)= [ (p{r,v,l -l')S{r,v,l')dl', 

Jo 

f°° 1 1 / s \ 

/ a(r,v,l) ^(r,v,t,l) dl = -C(f)[r,v,—V fi Cil}{r,v,s), 
Jo \ V \ V J 

or, equivalently, 

/ OO j-vt / ^ \ 

a{r,v,l) ^{r,v,t,l) dl = J — f2l,v,t - j dl. 

In the same way we obtain 


C 









ill, V, t 



dl, 


where the memory kernel is determined by Eq. (3.7) with p^ instead of p, while p^ is determined 
by Eqs. (3.6) with instead of a. 

As a result Eq. (3.2) is recast as the linear Boltzmann-like equation 


diJ-\-V ■ WJ = / K{r,v,v') 


— V J 4’{r,v,l)J^ — fil,v,t—-'^dl-\-F{r,v,t), (3.8) 


or concisely 

dtJ + v-VJ = iICM^-M)J + F, (3.8) 

where the scattering operator 1C is given by Eq. (2.3), and the memory operator is given by 


(j){r,v,l)'ip^r — fil,v,t -^ dl 

<P{r, V, t') J{r — vt', v,t — t') dt' (3.9) 

with 

v‘^(j){r,v,l) = F{r,v,t), 

I = vt, the operator JCi^ is determined by Eq. (3.9) with and instead of (p and F, respectively. 
The initial condition (2.4) results in the initial condition 

J\t=o = x{r,v). (3.10) 

Note that an equation similar to Eq. (3.8) (without absorption and the source term) was obtained 
in Ref. [20], see also [19]. 

If the extinction and scattering coefficients depend only on velocity, i. e., cr{r, v, 1) = (j{v) and 
(Tg(r, V, 1) = ajv), then the memory kernels are (p{r,v, 1) = (y{v) 5{l) and 4>s{r, v, 1) = crjv) 6{l). 
In this case Aiip = a{v)vj and = ajv)v'ip. Hence the Boltzmann-like equation (3.8) 

becomes the ordinary linear Boltzmann equation (without the acceleration term) 

+ V ■ VJ = / K{r,v,v')aJv')v''il>{r,v',t)dv' — a{v)v'ip + F{r,v,t). ( 3 . 11 ) 

Jy 

One can derive integral equations for the densities r] and J, see Appendix B. 



pvt 

MJ{r,v,t) = V 

Jo 
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4 Asymptotic solution for small mean free paths 


The Boltzmann-like equation (3.8) describes transport with an arbitrary distribution of the 
free path. But here we consider distributions with at least first and second finite moments. 

4.1 Assumptions 

We suppose that the “space” of velocities V is bounded and rotationally invariant, and 0 ^ V, 
i. e., all particles have nonzero velocities. We suppose also that the medium is isotropic and the 
kernel has the form 

K{r,v,v') = K{r, fi ■ fi',v,v') (4-1) 

where f2 = v/v, v = |u|, and the same for f2'. Besides, the kernel is supposed to be bounded 
from below and above: 

0 < Kmin ^ K ^ Ajnax < OO- 

This implies that jy Jy I A(-, u, v')\^ dv' du < oo. 

We suppose here that the scattering coefficient is equal to a^(r,v,l) = u!{r)a{r,v,l), where 
u is the scattering albedo. Then the Boltzmann-like equation (3.8) takes the form 

+ u • VV’ = (w/c - 1) Mil^ + F. (4.2) 

If the mean free path is small of order e, where e a small parameter, the extinction coefficient 
is represented as 

a{v,l) = , (4.3) 

see Appendix C. This implies representation for the memory kernel 

(j){v,l) = . (4.4) 

We suppose that the scattering albedo is given by 

a; = 1 — e^a(r). (4.5) 

This means that scattering dominates the transport, while absorption is small of order com¬ 
pared to scattering. Finally, we suppose that the source term is small of order e and represented 
by 

F = eF. (4.6) 


4.2 The equation and representation of the solution 

Taking into account all the assumptions we recast the Boltzmann-like equation as 

— fil, v,t — dl 

+ e'^F{r,v,t) (4.7) 

To find the asymptotic solution we represent the density in the form 

'ip{r, V, t) = V, t°) -I- 'ip^r, v, f) 

with t° = et and t' = t/s, where ■0° and are outer and inner solutions, respectively. The 
outer solution approximates the solution in the interior of the domain, while the inner solution 
approximates it in the initial layer. 


+ £v ■ V'ljj = [(1 — £^a)K, — X] 


rvt 


-</>{ V. 

eye 
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4.3 The outer expansion 

Substituting ■0° into Eq. (4.7) yields 




I 




el 


-(j)( V,- m r-f2l,v,t° - dl 


+ e F( r,v, — ] . (4.8) 


[This equation differs from Eq. (4.7) only in the first term on the left-hand side.] We assume 
that lias the asymptotic expansion 

V’° ~ V’o + + • • • 

Besides, we have the asymptotic expansion (see Appendix C) 


1 - 


I 


- 0 n, - ~ 




with 


(t>o ^ (t>o{v) = jf- and 4)^=(t)^{v) = 


(i) 2(1)^ 


^-1. 


(4.9) 


(4.10) 


Substituting the asymptotic expansions into Eq. (4.8) and equating coefficients of like powers of 
e yields the following equations: 


(fC-I) {(PqViIjq) = 0, 

{JC -1) {cPovri) = ^2■[vVro + {I- /Cl) (0in VV^g)], 
{JC - I) ( 0 OI 11 A 2 ) = d,.ro + v^2 ■ vri 


(4.11a) 

(4.11b) 


+ (X - /C) 


■ 'V'lpi + ^ % 


+ 1C (acpovipo) - F{r,v,— 


Eq. (4.11a) implies (see Appendix D) 


'4’oir,v,n = 


(4.11c) 


(4.12) 


where is a positive solution of the equation {I — IC)F = 0, and p is not yet determined. 

Remark: Instead of using the isotropic kernel (4.1) one may use a general kernel K{r, v, v'), 
impose a condition F = 1 and take (j)Q{v)v = Aq = const as in Ref. [10]. In this case 'ipQ{r,v,t) = 
p(r,t)/Ao. 

Eq. (4.11b) yields the solution (see Appendix D) 

V^?(r, n, t°) = -^[cF-n- V(Ap)] = ^ [cF - A^2 ■ V p - (f2 ■ VA) p ], 

cPqV (PqV 

where c = c{r,t°) is an arbitrary function, and 

A{r,v) = {X-lC,)-^{^jF^+^F. 

Substituting 'ipf into Eq. (4.11c) and using the solvability condition (see Appendix D) yields the 
equation 


Cd^op — div{B'Vp) — 'VB • Vp + I —AB + / aFv'^ dv \ p = — / F[r,v, — dn, (4.13) 


1 




47r 
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with 


f '1/v 

C{r) = -^dv 


and 


B{r) = \ [ -^^A{r,v)dv. 

37 

Hence i/jq is given by Eq. (4.12), where p satisfies Eq. (4.13). 

4.4 The inner expansion 

Since the outer solution satisfies the Boltzmann-like equation, the inner solution ijj^ should 
satisfy the homogeneous Boltzmann-like equation. Substituting ijj' into Eq. (4.7) with F = 0 
yields 


dpTp^ + ev ■ VV’^ = [(1 — e‘^a)JC — X] 


r-evt' 

Jo e 


I 


u, - V’M ^ -d/ 


ev 


(4.14) 


We assume that has the asymptotic expansion 

V'* ~ V’o + ^V’l + • • • 

Substituting this expansion and (4.9) into Eq. (4.14) and equating coefficients of e® yields the 
equation 

dpipl = {JC-I) [(PoViP'q) , 

which is in fact an ordinary differential equation. Its solution is 




(4.15) 


4.5 Initial conditions 

The initial distribution y does not depend on e, therefore we have the initial condition 

[i^o + i’o]\t=o = x{r,v). 

Due to properties of the scattering operator 1C (see Appendix D), in order to ensure that ipQ ^ 0 
as f —>■ 00 , the initial “value” c^qV '4>Q\^i_Q in Eq. (4.15) must be orthogonal to T in L 2 (V). This 

.,.0, 


results in the initial conditions 


it=o 


with 


and 


p\t=o = 


4>o{v)v 

ll^llL(v) 


V'o ^^^ = x{r,v) 


t=o 


(4.16) 

(4.17) 


4.6 Summary 

The asymptotic solution of the Boltzmann-like equation (4.2) with the extinction coeffi¬ 
cient (4.3) [this implies the memory kernel (4.4)], scattering albedo (4.5) and source term (4.6) 
is given by 

= 'ip^{r,v,t°) + 'ip}){r,v,f), 

with t° = et and f = t/e, where ^pQ is given by Eq. (4.12) with p obeying Eq. (4.13), and pjQ is 
given by Eq. (4.15). The initial condition (3.10) results in the initial conditions (4.16) and (4.17) 
for p and V’q, respectively. Note that ip^ decays exponentially as time increases. 
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5 The diffusion approximation to the one-speed Boltzmann-like 
equation 


5.1 The asymptotic solution 

The asymptotic solution, obtained in the previous section, becomes a diffusion approximation 
in the case of the one-speed equation. Since the velocity v is constant, the density can be given 

by 

Then the scattering kernel is given by 

K{r,f2,f2') = K{r,f2 ■ f2'), 
and the scattering operator becomes 


= [ K{r,f2 ■ 

The Boltzmann-like equation takes the form of the one-speed equation 


pvt 

Jo 


-|- vf2 ■ VV’ = (ojJC — X) 
where u) is the scattering albedo. The initial condition is 

^^ 14=0 = x{r,n). 


r — f2l, f2,t — — ] 


df2 + F{r, f2,t), (5.1) 


(5.2) 


Eq. (4.7) becomes 
edf-ip + evfi • V'i/’ 

= [{1 - e'^a)IC -1] 

We represent the density as 




dl 


-h e‘^F{r, 


i^ir, t) = f2, t°) -k f2, f) 


with t° = et and f = t/e, where V’° and ■i/’* are the outer and inner solutions, respectively. 
Eqs. (4.11) become 



{IC-I)ro = 0, 

(5.3a) 

e': 

cPov{JC-I)ri = n-[l + {l-K,)cP,]vVro, 

(5.3b) 


(PqV (IC - X) V’s = 9*0V’S + ■ VV'? + (Z - /C) (...) 



-k /C(a0owV’o) -F(r, f2, , 

(5.3c) 


with given by Eq. (4.10). The coefficient = K^{r) is the same as in Appendix D, but here 
it does not depend on v and v'. Eq. (5.3a) implies (see Appendix D, note that F = 1) 

ro{r,f2,n^roir,n, 

i. e., does not depend on f2. Eq. (5.3b) yields the solution (see Appendix D) 

ri{r,f2X) = c-Af2-Vro, 
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where c = c{r,t°) is an arbitrary function, and 


A{r) = — 




+ ’Pi 




(J") , K,{r) 


[ 2 ( 1 ) 


+ 


-2 ■ l-i^^(r.) 


Substituting V’l into Eq. (5.3c) and using the solvability condition yields the diffusion equation 


for ipQ 


with 


5^0V’S - div(E>° VV’S) + 


(5.4) 


Z)°(r) = 


v{J) \ (f) , K,{r) 


+ 


'2 ' l-K^{r) 


av 




3 [2(1) 

The initial condition (5.2) yields the initial condition 

^o\t=o = 

The term ipQ obeys the equation 

^hV’o = 

which implies the solution 


, K° =and /°(r,t°) = / Fir, 17,— dl7. 


(i) 




(5.5) 


The initial condition is 


V'o 


t=0 


x{r,^)-^[ x{r,^')df2'. 
47r 7§2 


(5.6) 


(5.7) 


In summary, the asymptotic solution of the one-speed Boltzmann-like equation (5.1), which 
satisfies the initial condition (5.2), is given by 

where ipQ satishes Eq. (5.4) and the initial condition (5.5), and is given by Eq. (5.6) with the 
initial condition (5.7). The term decays exponentially as time increases. 


5.2 The diffusion approximation 

Multiplying Eq. (5.4) by 47r and e yields the diffusion equation 

— div(F Vti) -|- Ku = /(r,t) (5.8) 

for the density of the particles 

u{r,t)= [ 'ip^{r,t°)dfl = 4TT'ip^{r,t°), 


where 


D{r) = 


{l)r 


{I") ^ K,{r) 


2{lf 


(1 — (jj)v , _ , 

K= - rn- and f{r,t) = 


(0 


F(r, 17, t) dl7. 


The initial condition (5.5) yields the initial condition 

[ X(r,f2)df2. 

J§2 


u\t=o = 


(5.9) 
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If fj = const, the distribution of the path length is exponential, i. e., p{l) = a and 
= a8{l). In this case the one-speed Boltzmann-like equation (5.1) becomes the ordinary 
one-speed linear Boltzmann equation 


+ vf2 ■ VV' = {(jjK, — Z) {avip) + F{r, 17, t), 
the moments are = 1/cr and = 2/cj^, and we obtain 


(5.10) 


(5.11) 


The diffusion approximation to the Boltzmann equation (5.10) can be obtained by the method 
of spherical harmonics [1, 3]. The latter yields the diffusion coefficient D = n/[3(T(l — 
which differs from the diffusion coefficient (5.11) by O(e^), since a; = 1 — O(e^). The absorption 
coefficient and initial condition, obtained by the method of spherical harmonics, coincide with 
the absorption coefficient (5.11) and the initial condition (5.9), respectively. 

Note that, if transport is considered in the two-dimensional space M^, the diffusion coefficient 


is 


D{r) 


2 l-K,{r) 


This coincides with the diffusion coefficient obtained in Ref. [27]: see Eq. (3) in Ref. [27] with 
p = 0 and IZ = K^. 
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A The phase space density, free path distribution and related 
functions 

The behaviour of the phase space density ^ along the characteristics of the hrst-order partial 
differential equation (2.1) can be completely described through the dependence on the variable 1. 
Therefore, if the dependence of ^ on position r, velocity v and time t is omitted this behaviour 
is described by the ordinary differential equation 

e' + ere = 0, / > 0, 

where ^ = ^{l) and a = (t{1) is the extinction coefficient. The solution of this equation is given 
by 

m = S{l)m, 

where ^ ^ 

5(/)=exp|—y a{l') = 1 — j p{l')dl'= J p{l')dl' 

is the survival probability, i. e., the probability that a free path of the particle is not less than I, 
and ^ 

p{l) = a{l) exp 1“ y ^i^') 
is the probability density function of a free path. 
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We define the memory kernel (p by [26] 




CS{X) l-Cp{X)' 


where 


r 

^/(A) = / 
Jo 

is the Laplace transform. This means that 


f{l)dl 


p{l)= [ cP{l-l')S{l')dl'. 

Jo 

If (T = const, i.e., the distribution is exponential, then (^(Z) = a6{l). This justifies the term 
“memory kernel”. 

The mean free path is given by 


roo roo 

{l)= / lp{l)dl= / 5(0 d/. 

Jo Jo 

The second moment is given by 

roo poo 

(0) = / l^p{l)dl = 2 lS{l)dl. 
Jo Jo 


B Integral equations for the densities 


Eqs. (2.2) and (3.5) (with the PDF p^ instead of p, corresponding to the scattering coefficient 
(Tg) yield the integral equation for the density p 


pir,v,t) = K, 


£‘ p.(..0 , (r - - I) di + p.(r... „«)x(r - „<.») 


+ F{r,v,t), 


cf. Eq. (2.8) in Ref. [5], where the initial distribution x is absent. 

Integration of the Boltzmann-like equation (3.8) along the characteristics and the initial 
condition (3.10) yield the integral equation 




= / {ICMs — M.) ip{r — vt',v,t — t')dt'+ / F{r — vt',v,t — t')dt'+ xif — 


C The asymptotic expansion of the memory kernel 


In this section the dependence of values and functions on v is omitted. 

We suppose that the mean free path (t^ is small of order e, where e a small parameter. To 
be specific, we represent it by (/) = sQJ), where (l) = 0(1) as e —)■ 0 [the symbol O is used 
here in the strict sense and means “strictly of the order of”, i. e., / = 0{g) means / = 0{g) and 
/ 7 ^ o(( 7 )]. Then the extinction coefficient is represented as 


This implies 


a{l) 




5(0 = s(l) 


and 
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with 


cm = 


Cp{X) _ \Cp{\) 


CS{\) 1 - Cp{\)' 

The first and second moments are given by 


and 


roo 

(l)= lp{l)dl 
Jo 

-9 

(i)= 

Jo 


We suppose that = 0(1) as £ —>■ 0. 

Expansion of the Laplace transform Cp{eX) in terms of e yields 

1 /t2\ 


Cp{eX) = 1 — (Z)eA + -(Z )e^A^ — ... 


Hence 

with 


C<p(^6X') — (pQ (p-^sX T (1)2^^X^ T • • • 


4’o — 7 t\ ‘Z^i ~ 




(l) 


- 1 


(an expression for 02 is not needed here). Therefore, we have formally the asymptotic expansion 

e V^/ 

(we consider the expansion in the weak sense). 


0( - ] ~ 0o'^(O + + 


D The Fredholm integral equation 

In this section the dependence of values and functions on r and t is omitted. 

Consider the Fredholm integral equation of the second kind 

{I-IC)C{v) = f{v), (D.l) 


where I is the identity operator, and 

}CC{v)= [ K{fl ■ fl',v,v')C{v')dv' 

Jy 

is the scattering operator with the kernel K(v,v') = K{f2 ■ f2',v,v'). The condition 

\K{v,v')\^ dv’ dv < oo 


means that 1C and its adjoint tC* are HilbertUSchmidt operators, which implies that they are 
compact in L 2 (V). The normalization 


[ K{v,v')dv = l ^ (I-JC*)1 = 0 

Jy 


(D.2) 


implies that C*('^) = 1 is the eigenfunction of the adjoint operator }C* corresponding to the 
eigenvalue 1. Therefore, the Fredholm alternative imposes the solvability condition 


f{v) dr? = 0. 


(D.3) 
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The strict positivity of the kernel K and Eq. (D.2) imply that the eigenvalue 1 is the largest 
in magnitude and simple eigenvalue of the operator JC, and the corresponding eigenfunction is 
positive, i. e., there exist a function > 0 such that 

{I-1C)^ = 0, (D.4) 

see Ref. [28]. Eigenfunctions, corresponding to other eigenvalues, are not nonnegative [28]. Due 
to isotropic scattering the eigenfunction ^ does not depend on the direction i7, i. e., ^(v) = ^(v). 
We need to solve the integral equation with 


f{v) = a(v) ■ n. 


(D.5) 


The solvability condition (D.3) is satished. To obtain a solution C, is expanded in the spherical 
harmonics 

oo n 

cw = E E Cniv)yri^) 

n=0m=—n 

with 

1 

CT{v)Yr{n)^Ci{v)-n 

m=—l 

(the functions are not the components of the vector-function The kernel is also expanded 
in the spherical harmonics 

OO n OO Oil 

K{o-o',v,v') = j2Kn{v,v') 

n=0 m=—n n=0 


where P„ are the Legendre polynomials. Due to orthogonality of the spherical harmonics we 
have 


/cic(v)yn^(^)] = i/c^c(v)]Y:^(f2), 


where 

= J K^{v,v')g{v'){v'fdv', 

and, therefore. 


(X-/C) [CAv)■fl] = {I-IC,)Cl{v)■^2. 


The operator (X — is dehned. Indeed, if (I — ICi) Ci(^) = 0, then Ci(^) ’ ^ = const P(r;). 

Therefore, Ci(^) = 0- 

As a result, we obtain that the solution of the integral equation (D.l) with the free term (D.5) 


IS 


({v) = ^^{v) + {X — JCi) ^ a{v) 




where c is an arbitrary constant. 
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